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Abstract. We extend our study [1] of the pion scalar radius in two-flavour lattice QCD to include two 
additional lattice spacings as well as lighter pion masses, enabling us to perform a combined chiral and 
continuum extrapolation. We find discretisation artefacts to be small for the radius, and confirm the 
importance of the disconnected diagrams in reproducing the correct chiral behaviour. Our final result for 
the scalar radius of the pion at the physical point is = 0.600 ± 0.052 fm^, corresponding to a value 

of (.4 = 4.54 ± 0.30 for the low-energy constant G of NLO chiral perturbation theory. 


1 Introduction 


Hadronic form factors carry crucial information about the 
internal structure of hadrons as bound states of quarks 
and gluons in Quantum Ghromodynamics (QGD). Since 
the dynamics underlying this structure is deeply non-per- 
turbative, the only known first-principles approach that 
allows for the extraction of form factors from QGD with 
full control of systematic errors is the use of lattice simu¬ 
lations. 

The scalar form factor of the pion, defined as 

(*5^) - (Pf)] '^ddd + m^uu |7r+ (pi)) (1) 

with the four-momentum transfer 


Q^ = -q^ = -{Pf-p^)\ (2) 


is not directly accessible to experiment for lack of a suit¬ 
able low-energy probe. However, the associated scalar ra¬ 
dius 
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can be related to the experimentally measurable cross sec¬ 
tion for TTTT scattering [5] using chiral perturbation theory 
(yPT). A notable feature of the scalar radius is that in 
yPT at NLO it depends only on a single low-energy con¬ 
stant, £ 4 , through 


^TT.phys 


( 4 ) 

where the pion decay constant is F = 92.2 MeV [5]. 

In lattice QGD, the scalar form factor of the pion is de¬ 
rived from a three-point function which receives contribu¬ 
tions from both quark-connected and quark-disconnected 
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diagrams (cf. figure . Due to the large numerical cost 
disconnected diagrams have often been neglected in lat¬ 
tice calculations. However, arguments based on partially 
quenched yPT indicate that quark-disconnected con¬ 
tributions to the pion scalar radius could be sizeable. Our 
earlier work [T], using lattice data at a single value of the 
lattice spacing and relatively large pion masses, has largely 
confirmed the prediction of yPT. In the present work we 
extend our study significantly, by including two additional 
lattice spacings, as well as lower pion masses. This allows 
for a combined chiral and continuum extrapolation of our 
results. In this way we obtain improved estimates for the 
pion scalar radius and the low-energy constant £ 4 , includ¬ 
ing a full assessment of systematic errors. 


2 Methods and Results 

2.1 Ensembles used 

In our calculation of the scalar radius we use Nf = 2 
dynamical flavours of non-perturbatively 0(a)-improved 
Wilson fermions. The gauge ensembles have been gener¬ 
ated by the GLS initiative and are listed in table 

Gompared to our previous work [T] we have included 
two more lattice spacings, a coarser one, a = 0.079 fm 
(labels A and B), and a finer one, a = 0.050 fm (labels N 
and O). For the intermediate lattice spacing, a = 0.063 fm, 
the G8 ensemble with a pion mass of 193 MeV has been 
included in our calculation. The use of different lattice 
spacings allows for performing a continuum extrapolation 
of our results for the scalar radius besides the chiral ex¬ 
trapolation. In general one also has to extrapolate to in¬ 
finite volume, however, all our ensembles fulfil m,rT > 4 
and thus we expect finite volume effects to be negligible. 
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Fig. 1. Quark flow diagrams for the connected and disconnected contributions to the scalar form factor of the pion. The cross 
represents the vacuum expectation value (tpip), which for the Wilson fermion discretisation used here contains a power-law 
divergent additive renormalisation. 


Label 

a[fm] 

lattice 

m.,r [MeV] 

m^L 

Ncig 

A3 

0.079 

96 X 48^ 

473 

6.1 

133 

A4 

0.079 

96 X 48® 

363 

4.7 

200 

A5 

0.079 

96 X 48® 

312 

4.0 

250 

B6 

0.079 

96 X 48® 

267 

5.1 

159 

E3 

0.063 

64 X 32® 

650 

6.6 

156 

E4 

0.063 

64 X 32® 

605 

6.2 

162 

E5 

0.063 

64 X 32® 

456 

4.7 

1000 

F6 

0.063 

96 X 48® 

325 

5.0 

300 

F7 

0.063 

96 X 48® 

277 

4.3 

351 

G8 

0.063 

128 X 64® 

193 

4.0 

348 

N5 

0.050 

96 X 48® 

430 

5.2 

477 

N6 

0.050 

96 X 48® 

332 

4.1 

946 

07 

0.050 

128 X 64® 

261 

4.2 

490 


Table 1. Overview of the CLS ensembles that have been used 
for the calculation of the scalar form factor of the pion. The 
lattice spacing given was determined using the 17 baryon mass 
[7]. Note that all ensembles fulfil > 4. 


2.2 Form factor calculation 


We explicitly calculate the connected and the disconnected 
contribution to the matrix element (7r| qq |7r) for three dif¬ 
ferent momentum transfers and determine the scalar pion 
form factor. The three momentum transfers are chosen 
such that the pion has momentum |pi| = 0 at the source 
and momentum |py| = 0, |p/| = or |pj| = 
at the sink. 

The quark loop required for the disconnected part is 
estimated using three stochastic sources per timeslice t 
and a generalised Hopping Parameter Expansion to 6th 
order. Further details can be found in [1]. 

We use ratios [5] of three- and two-point functions to 
extract the desired scalar matrix element from our data. 
For the connected contribution we use 


For the disconnected contribution we use smeared-smeared 
two-point correlation functions and the ratio 




C3{t,ts,Pz,Pf) 

C2{ts,Pf) 
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IC2{ts,Pf)C2{t,Pf)C2{{ts - t),Pi) 
C2its,Pi)C2{t,Pi)C2{{ts - t),Pf ) ' 


To increase the statistics for the disconnected contribution 
we average over four different pion source positions. 

In our previous work [1] we have discussed the re¬ 
maining time-dependences of these ratios due to periodic 
boundary conditions and the backward propagating pion 
in the two-point function. We take these dependences into 
account by dividing the ratios i?i and R 3 by the appropri¬ 
ate time-dependent factors. Furthermore, we have shown 
significant excited state contributions for small source-sink 
separations of ts < 1.5 fm. This has been confirmed on our 
extended set of ensembles. Thus, for all ensembles we only 
use data with ts > 1.5 fm. Further details on the extrac¬ 
tion of the scalar form factor can be found in ^ . 

For all ensembles we observe a significant contribu¬ 
tion of the disconnected diagram to the scalar form fac¬ 
tor. However, we find the relative disconnected contribu¬ 
tion to the form factor to be strongly dependent on the 
lattice spacing. Figure shows the disconnected contri¬ 
bution at vanishing momentum transfer divided by the 
corresponding connected contribution Fg’^(0)disc/^"s’^(0)con 
plotted against the pion mass. Different colours denote 
different values of the lattice spacing. Except for the light¬ 
est ensemble (G8), which might simply be an outlier, we 
find no significant dependence on the pion mass. How¬ 
ever, the relative contribution of the disconnected diagram 
has a very pronounced dependence on the lattice spacing. 
The smaller the lattice spacing, the smaller is the relative 
disconnected contribution {0 )disc/{0 )con to the form 
factor. 


Rl{t,ts,Pz,Pf) 


IC3{t,ts,p^,Pf)C3{t,ts,Pf,P^) 

C2{ts,Pi)C2{ts,Pf) 


(5) 

where the pion sink is located on ts and the scalar operator 
is inserted at time t. The pion source is placed at t = 
0 for simplicity. The required two-point and connected 
three-point functions have been calculated with Gaussian 
smearing p fTOlITT] applied to the source. 


2.3 The scalar radius 

The scalar radius of the pion is determined from the Q^- 
dependence of the form factor at vanishing momentum 
transfer. To estimate the scalar radius we use a linear 
parametrisation, i.e. 

F^{Q^) = F^{0) (^1 - i ^ ^ (7) 
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Fig. 2. The disconnected contribution F^(0)disc at = 0 
divided by the corresponding connected contribution. Blue, red 
and green points denote different values of the lattice spacing 
a = 0.050 fm, 0.063 fm and 0.079 fm, respectively. 


To determine the scalar radius on a given ensemble we 
fit a function of the form Q to our form factor results 
at the three different momentum transfers. Note, that the 
absolute normalization of the form factor is not known, 
since the multiplicative renormalisation Zg has not been 
determined. However, in the calculation of the scalar ra¬ 
dius Zg drops out. For the lightest ensemble G 8 , the ratios 
for the third momentum transfer were too noisy to obtain 
a useful signal. Thus, we have to resort to match a linear 
function to two values only for this ensemble. For all 
other ensembles where data for three different momentum 
transfers is available we find that the differences between 
using two or three momenta is not significant. 

To investigate a possible dependence of the scalar ra¬ 
dius on the lattice spacing we compare the results from 
ensembles with different lattice spacings at roughly the 
same pion mass. Within the CLS ensembles we have two 
sets of such ensembles (cf. table [^: « 270 MeV (B 6 , 

F7, 07) and « 325 MeV (A5, F 6 , N 6 ). In figure 
our results for the scalar radius on these ensembles are 
plotted against the squared lattice spacing a}. The upper 
panel shows results for a pion mass of TOtt ~ 270 MeV the 
lower panel for m.^. « 325 MeV. 

While for the smaller pion mass one can see a trend in 
the scalar radius for increasing lattice spacing, we do not 
observe such a trend for = 325 MeV. Although we do 
not find a significant dependence on the lattice spacing, we 
will include a term linear in c? for performing a combined 
extrapolation of the scalar radius to the physical point. 


3 Chiral and continuum extrapolation 

The use of gauge ensembles with different pion masses and 
lattice spacings allows for an extrapolation to the physical 
point vcIt^ —)■ m,r,phys and a —>■ 0. For the dependence on 
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Fig. 3. The scalar radius at fixed pion mass plotted against 
the squared lattice spacing a^. The upper panel shows results 
for a pion mass of « 270 MeV, the lower panel for « 
325 MeV. 


the pion mass we use the expression Q from NLO yPT. 
Since we use an 0(a)-improved action and because no 
operators that mix with qq at 0 {a) exist, we expect lattice 
artefacts to be of order a^. Thus, we include a term oc in 
a combined fit of the scalar radius for all gauge ensembles 
(cf. table [^: 


= - 
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+ba^ 


^ ( 8 ) 
The fit function ® has two ht parameters, the low-energy 
constant £4 and the coefficient b in the term which depends 
on the lattice spacing. When fitting our data for the scalar 
radius to the NLO expression Q we apply a cut in the 
pion mass by imposing < 500 MeV which excludes the 
ensembles E3 and E4. 

Additionally, we have repeated the fit using all data 
points (i.e. no mass cut) and a more aggressive mass cut 
m,r < 335 MeV. We find that the fit result is quite insen¬ 
sitive to mass cuts, and the fits produce consistent results 
with an acceptable x^/dof < 2. The spread among the 
results when different mass cuts are applied is negligible 
compared to the statistical error. Thus we use the fit with 
m,r < 500 MeV to quote a final result for the scalar radius. 

The result of this combined fit to our data is shown 
in figure where the scalar radius is plotted against the 
squared pion mass. The green, red and blue curves de¬ 
scribe the pion mass dependence at individual values of 
the lattice spacing, i.e. a = 0.079 fm, a = 0.063 fm and 
0.050 fm. Clearly, the lattice spacing dependence is very 
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Fig. 4. Our results for the scalar radius plotted against the 
pion mass including the result from a combined chiral and 
continuum extrapolation. Different colours denote different lat¬ 
tice spacings. The black line is the pion mass dependence in the 
continuum from our combined fit. The two right most points 
(open symbols) at « 0.37 GeV^ (E4) and m\ ~ 0.42 GeV^ 
(E3) are not included in the fit. The orange point at the phys¬ 
ical pion mass shows the value from TiTr-scattering [2]. 

small; we find b = 4.7 ± 12.5 in our fit, indicating that 
a possible dependence on cannot be resolved with the 
current statistical accuracy of our data. 

The black line in figure shows the pion mass depen¬ 
dence in the continuum, a = 0. At the physical pion mass 
we find for the scalar radius 

(r^)^ = 0.600 ±0.052 fm^. (9) 

This is in excellent agreement with the value 0.61±0.04 fm^ 
which was extracted from TTTr-scattering in [5] (cf. the or¬ 
ange point in figure]^. 

One has to note that the pion mass dependence of the 
scalar radius is completely driven by yPT (cf. equation 
W), and the fit parameter Z 4 only changes the offset of 
the black curve in figure We find that our results for 
the radius show the expected pion mass dependence from 
NLO yPT, confirming its applicability. 

From the combined fit we can extract the low-energy 
constant I 4 , which served as a fit parameter. We find 

I 4 = 4.54 ± 0.30 , (10) 

which is in very good agreement with the value I 4 = 4.62± 
0.22 quoted in the current FLAG report [12] for the Nj = 
2 flavour theory. 

In our previous work we found that the disconnected 
contribution to the scalar radius is important, which is 
in qualitative agreement with what has been found in 
partially quenched chiral perturbation theory |^. We can 
confirm this using our new data by repeating the same 
combined fit for the scalar radius using the connected 
contribution only: we indeed find a smaller value for the 


scalar radius = 0.532 ±0.042 fm^, correspond- 

mg to a smaller value for the low-energy constant = 
4.15 ± 0.24, which would disagree with the determination 
of £4 from other processes. 


4 Summary 

We have extended our previous study of the pion scalar 
radius using additional lattice spacings and pion masses in 
order to enable a fully controlled chiral and continuum ex¬ 
trapolation. We find that discretisation artefacts are mild, 
and that the pion mass dependence of the scalar radius 
is well described by NLO yPT. From a combined extra¬ 
polation to the physical point, we are able to extract the 
low-energy constant £4 = 4.54 ± 0.30 in good agreement 
with the FLAG [T2| average. 

The disconnected part of the pion three-point function 
contributes significantly to the slope of the scalar form fac¬ 
tor near = 0 and thus to the scalar radius of the pion. 
The inclusion of the disconnected diagrams is therefore 
essential in order to capture the correct physics. 
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